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Abstract — This paper addresses the problem of identifying 
the topology of an unknown, weighted, directed network run- 
ning a consensus dynamics. We propose a methodology to 
reconstruct the network topology from the dynamic response 
when the system is stimulated by a wide-sense stationary noise 
of unknown power spectral density. The method is based on a 
node-knockout, or grounding, procedure wherein the grounded 
node broadcasts zero without being eliminated from the net- 
work. In this direction, we measure the empirical cross-power 
spectral densities of the outputs between every pair of nodes 
for both grounded and ungrounded consensus to reconstruct 
the unknown topology of the network. We also establish that 
in the special cases of undirected and purely unidirectional 
networks, the reconstruction does not need grounding. Finally, 
we extend our results to the case of a directed network assuming 
a general dynamics, and prove that the developed method can 
detect edges and their direction. 

I. Introduction 

In recent years, complex dynamical networks have at- 
tracted considerable attention [1]. The power grid, the In- 
ternet, the World Wide Web, as well as many other bi- 
ological, social and economic networks [2], are examples 
of networked dynamic systems that motivate this interest. 
The availability of datasets describing the structure of many 
real-world networks has allowed to detect the presence of 
common patterns in a large variety of networks [3], [4]. 
In this paper, we address the problem of reconstructing the 
structure of an unknown network of dynamical nodes from 
observations of its input-output behavior. 

The problem of network reconstruction is crucial in a wide 
variety of disciplines such as biology [5]-[8], physics [9]- 
[11] and finance [12]. A wide collection of approaches have 
been proposed to solve the network reconstruction problem. 
For example, we find in the literature several papers that 
approach this problem using an optimization framework, 
such as [8], [11], [13]. In these papers, the reconstruction 
problem is stated as the optimization problem of finding 
the network that maximizes a function that measures the 
sparsity of the network (e.g. ^i-norm) while conforming 
to known a priori structural information. Although the as- 
sumption of sparsity is well justified in some applications 
(e.g. biological networks), this assumptions might lead to 
unsuccessful topology inference in other cases, as illustrated 
in [14] and [15]. When the unknown network is known to 
be a tree, several techniques for network reconstruction were 
proposed in [12], [16] and [17]. More recently, Materassi and 
Salapaka proposed in [18] a methodology for reconstruction 

^The authors are with the Department of Electrical and Systems Engi- 
neering, University of Pennsylvania, Philadelphia, PA 19104-6228 USA. 
{shahin, preciado}@ seas . upenn . edu 



of directed networks using Wiener filters. Although effective 
for many networks, this methodology is not exact when 
two nonadjacent nodes are connected to a common node 
with directed edges pointing towards the common node. In 
[19], Nabi-Abdolyousefi and Mesbahi proposed a technique 
to extract structural information, such as node degrees, 
of an undirected network running a consensus dynamics. 
Subsequently, they find a collection of undirected graphs that 
are consistent with this structural information. Furthermore, 
the work in [20] provides a method that performs a complete 
identification of undirected networks via a procedure called 
node knock-out. 

In this paper, we propose an approach to reconstruct 
the structure and weights of a directed network from the 
output of an agreement dynamics run on the network. We 
also revisit the problem of network reconstruction for more 
general dynamical systems. We develop methodologies to 
unveil the network structure from the dynamic response of 
the network when the system is driven by stochastic inputs. 
More specifically, we assume that the system is stimulated 
by a collection of wide-sense stationary noises with unknown 
power spectral densities. Considering several cases, we pro- 
pose methodologies to recover the network topology from 
the empirical cross-power spectral densities of the outputs 
between every pair of nodes. We first consider the case 
of undirected networks running a consensus dynamics, and 
propose an algorithm to reconstruct its unknown topology. 
Using the node knockout procedure proposed in [20], we 
extend our results to the directed networks. In this scenario, 
the node knockout is equivalent to state grounding, where the 
node broadcasts a zero state to its neighbors without being 
eliminated from the network. We also prove that for purely 
unidirectional networks (networks with no reciprocity [21]), 
there is no need to run the grounded consensus to perform the 
network identification. Finally, we consider the reconstruc- 
tion of directed networks assuming a general dynamics. We 
establish that, without the knowledge of the power spectral 
densities of the input vector, a Boolean reconstruction is still 
possible, i.e., we can detect edges and their direction, but not 
their weights. 

The rest of the paper is organized as follows. In section 
II, we provide some nomenclature needed in our exposi- 
tion and describe the network reconstruction problem under 
consideration. Assuming the system is driven by wide-sense 
stationary noise, in section III. A, we provide a technique 
to recover the structure of an undirected network running 
consensus dynamics. Section III.B covers the extension of 
the problem to a directed network. Section III.C addresses 
the reconstruction of directed networks following general 



dynamics. Section IV concludes. 

II. Preliminaries & Problem Description 

In this section, we introduce a series of definitions used 
throughout the paper. Let Q = (V, £) be an unweighted, 
undirected graph, where V = {i>i, . . . ,v n } denotes a set of 
n nodes and £ C V x V denotes a set of in undirected 
edges. If {vi,Vj} 6 5, we call nodes Vi and «j adjacent (or 
first-neighbors), which we denote by V{ ~ w,-. A weighted, 
undirected graph is defined as the triad W = (V, £ , J 7 ), 
where V and £ are the sets of nodes and edges in W, and 
the function T : £ — >M associates real weights to the edges. 
Similarly, a weighted, directed graph is defined as the triad 
T) = (V, £d, J-d), where V is the set of nodes and £d is the set 
of directed edges in V, where a directed edge from node Wj 
to node Vj is defined as an ordered pair (Uj, vj). Furthermore, 
Td is a weight function £d ~> 

In an unweighted, undirected graph Q, the degree of a ver- 
tex Uj, denoted by dcg (vi), is the number of nodes adjacent 
to it, i.e., dcg(wj) = G V : {^Vj} G £}\. This defini- 
tion can be generalized to both weighted and directed graphs. 
For weighted graphs, the weighted degree of node is equal 

to deg (v^ = J2 j:{vi.Vj}e£ ^~ ({ Vi ' V j})> ^' e -' SUm °f '' rie 

weights associated to edges connected to Wj. For weighted, 
directed networks, we define the weighted in-degree of node 

Vi as deg in = E,j;( Vj , Vi )es d (i v j> v i))- 

The adjacency matrix of an unweighted, undirected graph 
Q, denoted by Ag = [ffly], is a n x n Boolean symmetric 
matrix defined entry-wise as Oy = 1 if nodes Vi and Uj are 
adjacent, and = otherwise. We define the Laplacian 
matrix Lg of a graph Q as Lg = Dg — Ag where Dg is 
the diagonal matrix of degrees, Dg = diag ((deg (i>i))iLi)- 
For simple graphs, Lg is a symmetric positive semidefinite 
matrix, which we denote by Lg > [22]. Thus, Lg has 
a full set of n real and orthogonal eigenvectors with real 
nonnegative eigenvalues = Ai < A2 < ... < A„. 

Similarly, the weighted adjacency of a weighted graph W 
is defined as Aw = [iVij], where Wij = J- ({vi,Vj}) for 
{vi,Vj} G £, and Wij = if {w^fj} ^ £. We define 
the degree matrix of a weighted graph W as the diagonal 
matrix Dyy = diag ((deg (vi))™ =1 ). The Laplacian matrix 
of a weighted, undirected graph W, is defined as Lyv = 
Dy\; — Ay^. Furthermore, the adjacency matrix of a weighted, 
directed graph V is defined as A-p = [dy], where dy = 
Fd{{vj,Vi)) for G and dy = if (vj,Vi) £ 

£d- We define the in-degree matrix of a directed graph T> 
as the diagonal matrix Dx> = diag ((deg in (vi))2=i)- ^he 
Laplacian matrix of I? is then defined as L-p = Dx> — A%>. 
Note that the Laplacian matrix satisfies Lgl — Lyyl = 
Lp\ = 0, i.e., the vector of all ones is an eigenvector of the 
Laplacian matrix with corresponding eigenvalue 0. 

In this paper, we focus on the continuous-time non- 
autonomous model of consensus dynamics in directed net- 
works, described as 



where Lx> is the Laplacian matrix of a weighted, directed 
network V; and the vectors x(t),w(t) ,y(t) G M" are the 
state, input, and output vectors, respectively. We assume that 
we do not have access to the network structure, i.e., L-p 
is unknown. In this context, we consider the problem of 
identifying the topology of the directed network V when 
the (cross-)power spectral densities of the output vector 
y(t) are measured empirically while the stochastic input 
vector w (i), injected at the nodes of the network, has 
unknown power spectral characteristics. We also investigate 
the reconstruction problem in the case that ([TJ follows a 
general dynamics Gp in lieu of — Lp. In any case, we 
assume the input w (t) = [wt (t)] is a vector of uncorrected 
wide-sense stationary processes. 

Definition 1 (Wide-Sense Stationary): A continuous-time 
scalar random process w(t) is wide-sense stationary (WSS), 
if it satisfies the following properties: 

PI. (j, w (t) = E(tu(i)) = fi w (t + t) for any r G K. 

P2. R w (h,t 2 ) 4 E(w(ti)wfo)) = R w {h +T,ta + 
t) = R w (h - t 2 , 0) for any rel. 

The reconstruction methods proposed in this paper take 
the output vector y (t) = [iji (t)] and deliver the network 
structure as a function of the (cross-)power spectral densities. 

Definition 2 ((Cross-)Power Spectral Density): The 
cross-power spectral density of two WSS signals, t/j (t) and 
yj (t), is the Fourier transform of their cross-correlation 
function, i.e., 

W<") =HRy iyi (r) =E( yi (t) yj (t-r))}, 

where T {■} is the Fourier transform operator. The power 
spectral density of yi (t) is defined as 

S V M ± T {Ryt (r) ^ E( yi (t) yi (t - r))}. 

In order to detect the links through observation of outputs, 
we need to ascertain that system dTJ is driven by a collection 
of noises with nonzero power spectral densities. To capture 
this idea, we commence with the following definition: 

Definition 3 (Excitation Frequency Interval): The 
excitation frequency interval of a vector w (t) of wide-sense 
stationary processes is defined as an interval (—0, fi), 
with £1 > 0, such that the spectral densities of the input 
components t0{ (t) satisfy S Wi (uj) > for all ui G (—CI, f2), 
and all i G {1, 2. n}. 

We end this section by stating our assumptions. Through- 
out the paper we impose the following conditions on the 
input vector: 

Al. The collection of signals {wi(t),i = l,...,n} are 
uncorrected zero-mean WSS processes such that, 
for any t, r € M, 

E(t« i (t)«> i (T))=0, for 1^ j, 

and 



x(t) = -L v x{t) + w (i) , y{t) = x(t), (1) 



R Wi (T)=E{w i (t)w i (t + T))=R w (T). 



A2. There exists a nonempty excitation frequency inter- 
val 

III. Main Results 

We now consider several cases and present methodologies 
to reconstruct the structure of an unknown network from 
observations of its temporal response ([T}. We provide re- 
construction techniques for weighted, undirected networks 
(Sect. III. A), and weighted, directed networks (Sect. III.B) 
following consensus dynamics. We briefly, revisit the recon- 
struction of directed networks following general dynamics 

(Sect. m.c). 

For the consensus dynamics (as well as general dynam- 
ics) of weighted, directed networks, the following lemma 
provides an explicit relationship between the cross-power 
spectral densities of two outputs, y^ (i) and yj (t), and the 
power spectral density of the input w k (i) (which we assume 
to be identical, i.e., S Wk (uj) — S w (uj) for all k): 

Lemma 4: Given assumptions (A1)-(A2), the following 
identity holds 



= ef (uj 2 I - ju:(L v - L v ) + L l v L v )- L e^ (2) 



S w (u) 

for any ui 6 (— 0, fi), where 

S ViVj( u ) = F{ R mDi( T )} and S w (uj) = F{R w (t)}, 

for any 1 < i,j < n. 

Proof: The transfer function corresponding to the state- 
space equations ([]]) is H (uS) — (juil + Lp)^ 1 . The transfer 
function from the k-th input w k (t) to the i-th output yt (t) 
is defined as 

H ki (uj) = efijujl + Lp)~ 1 e k . 

Therefore, the power spectral density of the i-th output j/i (t) 
when the input is w (t) = w k (t) e k (i.e., a WSS noise on 
the fc-th node) is equal to 



Syi(u) = Hki(u))HZi(u)S Wk (ui). 



(3) 



On the other hand, the transfer functions from input 
w k (t) to the outputs yi (t) and yj (t) are, respectively, 

Y, (w) /W k (w) = H kt (uj) and Yj (w) /W k (uj) = H kj (u). 
Hence, Yj (uj) jYi (uj) = H^(uj)H k j(uj) which implies 



(4) 



Since S Wk (w) — S w (u>) for all k, we can combine (f3]) and 
(|4]l to obtain 



S ViV 3 { UJ ) = H kl (uj)H kj (uj)S w (uj). 



(5) 



When the input is w (t) = Y^ k =i w k (t) &k> with 
¥,(wi(t)wj(T)) = for i j, we can superpose (0 over 
1 < k < n, to obtain the following for any ui G (— f2,fi): 



n 

= ^ ef (-jwJ + L T >y 1 e k e[ Qui + L B ) _1 e fc 
fc=l 

n 

= ^eJ(-jcjJ + L-D) _1 e fc efc(jaj/ + L5) _1 e l 
fc=l 

n 

= e J(-ja;J + L-D) _1 (^e fc efe)(jcj/ + L^)" 1 ^. 

fc=i 

As X)fc=i e fc e fc = we can simplify the last equation to 

= ef (uj 2 I - juLv +ju;Lv + L^,L v y 1 e J , 

which is the desired relation. ■ 
Corollary 5: Given assumptions (A1)-(A2), and substitut- 
ing — Lp in (UJ with any negative semi-definite G-p, the 
input-output power spectra relationship is as following 



S w (uj) 



ef (a; 2 / - fa(G% - Gp) + G^G^ej, 



for any ui G (— f2, 51). 

Proof: In the proof of Lemma |4] we did not use any 
properties of Laplacian, so we only need to replace — Lp by 
any negative semi-definite Gp in © to derive the result. ■ 
In our analysis in the next subsections, we frequently 
invoke the result of Lemma |4] and Corollary [5] 

A. Undirected Laplacian Identification from Stochastic In- 
puts 

In this subsection, we propose an approach to reconstruct 
the topology of an unknown weighted, undirected network 
W from the output of (|T). The formal statement of the 
reconstruction problem for weighted, undirected network can 
be stated as follows: 

Problem 6: Consider the dynamical network in ([T), with 
L-p = Lyv where W is an unknown weighted, undirected 
graph. Find the structure of W, from the empirical (cross- 
)power spectral densities^ of the outputs, i.e., S Vi (ui) and 
S Vi y .(ui) for all 1 <i,j < n. 

In the case of weighted, undirected networks, we can use 
the result of Lemma|4]in the form of the following corollary: 

Corollary 7: Consider the network dynamics (fl3, when 
V is a weighted, undirected network W. Given assumptions 
(A1)-(A2), the (cross-)power spectral densities satisfy 



J ViVi 



S w (u) 



ej {w 2 I + L 



w) 



(6) 



for any 1 < i,j < n. 

Proof: Since Lyv is symmetric, replacing L-p in (O by 
Lw, the proof follows immediately. ■ 

'There are several methods to empirically measure the (cross -)power 
spectral densities (see e.g., Welch's method [23]). 



We now proceed with a theorem that can be used to 
reconstruct a weighted, undirected Laplacian Lyv> from 
the empirical (cross-)power spectral densities S Vi (uj) and 
S Viy .(ui) for all 1 <i,j < n. 

Theorem 8: Consider the network dynamics ([TJ, when T> 
is a weighted, undirected network W. Let us define the 
matrix of cross-correlations as S (cj) = [S yiyj (io)\- Then, 
given assumptions (A1)-(A2), we can recover the weighted, 
undirected Laplacian Lyv as 



Ly\) — LO 



i 



1/2 



(7) 



for any i, and for any lo in the excitation frequency interval 
(-f2.fi). 

Proof: According to Corollary [7] we have 



S(u) = S w {u){u 2 I + L 2 w )-\ 



which yields 



Hence, 



S w (uj) 



(lu 2 I + L 2 w ). 



r 2 



SMS-^co)-^ 2 !. 



(8) 



(9) 



We can derive an expression for S w (lo) in terms of S 1 (oj) 
and lo, as follows. Post-multiplying © by the vector 1, we 
get 

, ,2 



1 



S w {u)) 



(u 2 I + L 2 w )l 



S w (io) 



since L 2 *,! — 0. Therefore, 



S w (ui) 



(10) 



(ID 



for any i. Substituting ( fTTT i in (O, we obtain the statement 
of our Theorem. (Notice that Lyy is the Laplacian of an 
undirected network, hence, it is positive semidefinite). ■ 
In practice, since [S _1 (u;)l] . is a noisy measurement, 
according to ( TTOb one can average out noise by replacing 
[S _1 (cj)l] with l T S _1 (w)l/n, to achieve improved em- 
pirical results. Hence, (0 can be rewritten as 



Lyy — LO 



1/2 



Remark 9: To derive node degrees computationally effi- 
cient, one can circumvent the implicit eigenvalue decompo- 
sition involved in the result of Theorem^] for computing Lyy 
from L^y, since 



[L'w\u — [{D w - A w ) 2 \u 



21 Vh 

[Dyy - D W A W - A W D W + A w ]i 



= deg 2 (vi) + deg(vi). 

Therefore, computing L 2 ^, we obtain the degrees by solving 
the quadratic equation above for any i. 



B. Directed Laplacian Identification from Stochastic Inputs 

In this section we address the problem of reconstructing a 
weighted, directed Laplacian Lp when the input is a vector 
of uncorrected wide-sense stationary processes w (t). 

Problem 10: Consider the dynamical network in ([TJ, 
where T> is an unknown weighted, directed graph. Recon- 
struct V, from the empirical (cross-)power spectral densities 
of the outputs, i.e., S Vi (Lo) and S yiyj (Lo) for all 1 < i,j < n. 

In what follows, we propose a reconstruction approach 
based on a grounded consensus dynamics, similar to the one 
proposed in [20] for the reconstruction of undirected graphs. 
The grounded consensus is defined as follows: 

Definition 11 (Grounded Consesus): The consensus dy- 
namics with node Vj grounded takes the form 



x(t) = -L-Djxit) +w(t) , y(t) = x(t), 



(12) 



where L-p. £ ]R(' l_1 ) >< ( Tl_1 ) is obtained by eliminating the 
j-th row and the j-th column from Lx>. 

The consensus dynamics with node Vj grounded describes 
the evolution of the network when we force the state of node 
Vj to be Xj(t) = 0. We now state a lemma and a corollary 
to extract LT^Ld from the ungrounded consensus ([TJ, and 
Z7p Lxij from the grounded consensus (fT2l , respectively. 

Lemma 12: Consider the network dynamics ([TJ, when T) 
is a weighted, directed network. Let us define the matrix 
of cross-correlations as S (a;) = \S yiVj {lo)\. Then, given 
assumptions (A1)-(A2), we can compute L^,Lx> and Lx> — 
Lf*, as 



L>t>Ltj 



LO' 



Ln 



RejS-V)} 
[RelS- 1 ^)}!], ~ 
ImjS-^q;)} 
[Re{S-i( W )}l] 



for any i, and for any lo in the excitation frequency interval 

Proof: Let us consider the matrix S(lo) = [S yiyj (Lo)}. 
According to (O, we can separate the real and imaginary 
parts of its inverse as 

Re{S-V)} = ^-^(w 2 / + L%L V ) (13) 

Im{S-\Lo)} = -—^-(LoL v -LoLl). (14) 

O w (LO) 

Post-multiplying ( [T3l by 1, and taking into consideration that 
L^L-pl — 0, we obtain 



S w (lo) = 



[Re{S-i(w)}l] i 



(15) 



for any i. Plugging (H~5b into ( fOl ) and ( TT4b . the proof follows 
immediately. ■ 
Corollary 13: Consider the network dynamics (fl~2), where 
Li) is the grounded Laplacian matrix of the weighted, 
directed network V grounded at vj. Let us define the matrix 
of cross-correlations as S (lo) — [S ViVk (Lo)} for all i,k ^ 



j. Then, given assumptions (A1)-(A2), we can compute 

LL\.L-p. as 

for any i ^ j, and for any to in the excitation frequency 
interval (—0,0). 

Proof: By Corollary [5] it holds that 



SyiVk ( w ) 



el(uj 2 I-juj(L Vj - L T V ) + Ll L Vj )- x e k 



S w (lo) 

Since S(cj) = [S ViVk (ui)], we have that 
2-1 



1 



Re{S-»} = — ^(uj 2 I + LlL Vi ). (16) 
S w [w) 3 

Substituting the expression for S w (w) from Lemma [TU we 
have 



L vM=« I [Re{S- 1 (w)}l] 4 ~ 



In general, the results of Lemma Q~2] might not be in- 
formative enough to extract the underlying structure of the 
network without running the grounded consensus. We will 
see at the end of this section how running both ungrounded 
and grounded consensus will lead to exact reconstruction 
of the network. However, imposing certain conditions on 
the network graph, allows us to perform the reconstruction 
without recourse to grounding. For instance, for a network 
that does not contain any bidirectional edge, i.e., if (vj,Vi) 6 
Ed, then Vj) ^ £d, we can employ Lemma [T2l to identify 
the network. The adjacency matrix of a purely unidirectional 
network satisfies 



tv(A 2 v ) = 0. 



(17) 



We now proceed to the identification technique of a purely 
unidirectional network. 

Corollary 14: Suppose the conditions in the Lemma Q~2] 
hold. Also, suppose the weighted, directed network satisfies 
dlTt . Then, the entries of adjacency are recovered as 



max < ujn 



( [ImjS- 1 ^)}]^ 
^lTRejS- 1 ^)}! 







Proof: For any k, by (fT~5T > we have 

[Re{S-V)}lU = -(^RejS- 1 ^)}!). 

n 

Therefore, under purview of Lemma [T2l we obtain 

/ ImjS- 1 ^)} 



Li 



Lv- wn ^!T Re{S -i (a;)}1 



If the ij-th entry of the right hand side in the above equation 
is negative, since there is no bidirectional edge, [Ax*]^ ^ 
and [At>] = 0. If the ij-th entry of the right hand side 
in the above equation is positive, the ji-th entry would be 
negative which implies [^Id]^ ^ and [^4©]^ = 0. In case 



the ij-th entry is zero, no directed edge between vi and Vj 
exists, and the proof is complete. ■ 

In the next theorem, we show using grounding, there is 
no need for structural conditions to solve Problem [TOl 

Theorem 15: Consider the network dynamics (Q3, where 
V is a weighted, directed network. Let us also consider 
the grounded consensus dynamics in (1121 1. when node Vj 
is grounded. Then, given assumptions (A1)-(A2), we can 
recover the entries of the weighted, directed adjacency matrix 
Ad as 



[Ma = 



[L>v L i 



L'L\ j L- Dj 



for i < j, 
for i > j, 



for any i ^ j. 

Proof: For simplicity, we consider the case j — n (for 
any j ^ n, we can transform the problem to the case j = 
n via a simple reordering of rows and columns). Running 
the ungrounded consensus (Q} and applying Lemma [12] we 
recover LI^Ld- Also, running the grounded consensus (TT~2T > 
and applying Corollary Qj] we can recover Lp n Ip n . Hence, 
for any i < n 

[LvLv^i-lLTp^Lv^ii = ^2 [-kp]fci~yi [ L v] ki = [Lvf ni ■ 
Hence, we can recover the adjacency as 



Lj, L Vj 



for any i < n. The same analysis holds for any other j ^ n. 

■ 

C. Boolean Reconstruction of General Dynamics 

In this section, we study the reconstruction of a general 
dynamics G-p in weighted, directed networks, characterized 

as 

x(t) = G v x{t) + w(i) , y(t) = x{t), (18) 

where Gx> is a negative semi-definite matrix. 

It turns out that for the general dynamics (fT8~t . complete 
identification including weights is not possible. However, we 
can propose an algorithm for a Boolean reconstruction of the 
system which detects the existence and direction of edges in 
the directed network ( fT8l . 

Problem 16: Consider the dynamical network in dTBT l. 
where T> is an unknown weighted, directed graph. Recon- 
struct an unweighted version of V, from the empirical (cross- 
)power spectral densities of the outputs, i.e., S yi (uj) and 
S ViVj {u}) for all 1 < i, j < n. 

To solve Problem [16] similar to Definition QT| we define 
the Grounded Dynamics at node Vj as 



b(t)=G Vj x(t)+w(t) , y(t)=x(t), 



(19) 



where Gx> G ]R( n_1 ) >< ( n_1 ) is obtained by eliminating the 
j-th row and the j-th column from Gp. We now state a 
theorem to reconstruct an unweighted version of T>. 



Theorem 1 7: Consider the network dynamics ([T8V where 
I? is a weighted, directed network. Let us also consider 
the grounded consensus dynamics in ( fT9] l, when node Vj 
is grounded. Then, given assumptions (A1)-(A2), we can 
recover a scaled version of off diagonal entires for matrix 
Gt> as 



^[Re{[S- 



-[s- 



i}] 



for i < j, 



S w [Re{[S^] u - [S- 1 ],-!,,-!}] for i > j, 



for any i ^ j. 

Proof: Without loss of generality, let j — n as in the 
proof of Theorem Q3] By Corollary |5] it holds that 

S^HRelS-V)} = uj 2 I + GlG v 
S w (u J )Re{S- 1 (uj)} =uj 2 I + GI G Vn , 



which entails, for any i < n 

[ G v)ni = [G%Gv]i 

= S w {lo)([Rz{S-\u)} 



p*e{S-»}]< 



which implies the result of the theorem. The proof for j ^ n 
follows with similar arguments. ■ 

Remark 18: Results of previous sections regarding Lapla- 
cian (Sect. III.A-B), could be extended to any negative semi- 
definite Gt> with a nonzero nullity, as long as a nonzero 
vector in the kernel of Gx> is known. 

Theorems [8] Q3] and [17] are functional for any frequency 
w ^ 0, and we do not require the knowledge of entire 
power spectral densities. In fact, we only need the power 
spectral densities evaluated at one frequency, which dramat- 
ically reduces the computational complexity. One can also 
evaluate the spectra at several frequencies to average out the 
measurement noise and make a robust reconstruction. 

IV. CONCLUSIONS 

In this paper, we addressed the problem of identifying 
the topology of an unknown, weighted, directed network 
running a consensus dynamics when the system is stimu- 
lated by a wide-sense stationary noise of unknown power 
spectral density. We proposed a methodology for network 
reconstruction based on sequentially grounding nodes in the 
system. We showed how to reconstruct the topology of the 
network from the empirical cross-power spectral densities of 
the outputs between every pair of nodes. We also established 
that, in the special cases of undirected (i.e., perfect edge 
reciprocity) or purely unidirectional networks (i.e., no edge 
reciprocity), our reconstruction procedure does not require 
any grounding. Finally, we extended our results to the case 
of a directed network with a general dynamics, and proved 
that the developed method can detect the presence of edges 
and their direction. 
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